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Abstract 

For a second order differential operator A(x) = — Va(x) V + b'(x) V + V(b"(x) ■) 
on a bounded domain D with the Dirichlet boundary conditions on dD there 
exists the inverse T(\, A) = (XI + A)' 1 in Li(D). If pi is a Radon (probability) 
measure on Borel algebra of subsets of D, then T(\,A)fi G L p (D),p 6 [f, d/(d — 
1)). We construct the numerical approximations to u = T(\,A)p in two steps. 
In the first one we construct grid-solutions u n and in the second step we embed 
grid-solutions into the linear space of hat functions u(n) £ Wp(D). The strong 
convergence to the original solutions u is established in L P (D) and the weak 
convergence in Wp(D). 

AMS subject classification: (2000) 35J20, 35J25, 35J15, 65N06, 65NI5 
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1 Introduction 

If the tensor- valued function {aij}ff (diffusion tensor) on a bounded domain Del 1 ' 
satisfies the strict ellipticity conditions, then the second order differential operator 
in divergence form: A — —J2ij di a ijdj + Y^i^'fli + with homogeneous 

Dirichlet boundary condition on dD, has an inverse A~^ in terms of an integral op- 
erator mapping the Radon measures on Borelian sets S(-D) into Wp(D) for each 
p 6 [1, dj (d — 1)) |BQ[ ILR2] . There are partial results on numerical solutions of the 
boundary value problem Apii = \x [ClllLR2[|LR3| . Here we extend results of |LR3j to 
a general boundary value problem on D C R. d . The operator Ad is discretized by a 
system matrix A n and numerical solutions are represented by grid- functions u„ . Then 
u„ is imbedded into the linear space of hat functions and compared with the solution 
of the original problem in order to prove convergence. Our particular intention is 
to construct numerical approximations u„ with system matrices A n possessing com- 
partmental structure, that is a matrix having positive diagonal elements, non-positive 
off-diagonal elements and non-negative column sums. Apparently, such matrices are 
transpose of M-matrices. Because of our determination to look for approximations A n 
with the compartmental structure, we have faced a number of non-typical problems in 
numerical analysis. All of them stem from the construction of the numerical schemes, 
rather then from the structure of convergence proofs. 

Various definitions and auxiliary results which are necessary to prove convergence 
are given in Sections [2] and [3l In the case of dimension d = 2 the proposed numerical 
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scheme is already on a level of an algorithm and can be straightforwardly applied. 
The construction of discretizations A n with the compartmental structure for any 
dimension is carried out in Section 21 Grid-solutions of discretized problems are 
imbedded in the linear space of hat functions and the obtained approximate solutions 
are analyzed from the standpoint of convergence. The convergence in W^-spaces is 
proved in Section[5]and the convergence in W^-spaces in Section[6j Section[7]provides 
two examples, demonstrating the efficiency of constructed schemes. 

2 Definition of the problem 

Elements of M. d are denoted by x, m, p etc. The Euclidean norm in IR d is denoted by 
| ■ | . The only open subsets of M. d considered in this work, are bounded and connected 
open sets with Lipshitz boundary [Ma, StcJ . We call a subset of this kind a domain 
with Lipshitz boundary. We denote it by D, and its boundary by 3D. For a set S C M. d 
the closure is denoted by S and sometimes by cls(S). 

The Banach spaces of functions (R d ), C (k *> (D) are defined as usually. Their 
norms are denoted by || • ||qo . The closure of functions in C^(M. d ) with compact 
supports determines the subspace Cg^(M d ). The closure of functions in C(D) with 
supports in D determines the subspace C(D), and then (7 (fe) = C (fe )(IX) n C(D). 
The Lp-spaces as well as Sobolev -spaces are defined in a standard way |Ma|IStej . 
Their norms are denoted by || • || p and || • || P) i, respectively. For each p, 1 < p < oo, 

the norm of W*(D) is defined by ||u|| p ,i = (\\u\\ 2 p + || Vu^) 1/2 , where || Vu || p = 
( Y^j=i II ®j u lip) ■ Because the domain D has Lipshitz boundary, the space (D) 
can be realized as the space of continuous functions on D, for which the first partial 
derivatives are elements of L oa (D). The completion in the norm of W p (D) of functions 
in CW(.D) n C(D) is denoted by W${D). The linear space C(D) n W${D) is dense 
in Wp(D) for all p, 1 < p < oo, and the closure of C{D) n Wp(D) in the norm of 
Wp(D) is equal to W£(D) [Ma]. Let X be a Banach space and its dual. Then 
the value of / G at u G X is denoted by (f\u). Let 11(D) be the convex set of 
positive Radon measures /i on Then (v \ n) = J D v(x)/i(dx) is well defined for 

v G W^ C (D), We say that a sequence of \x n G 1Z(D) converges weakly to /i S 1Z(D) if 
lim n ( v\fi n } = (v\n ) for each v € C(D). 

Let Is be the indicator of S C M. d , i.e. ILs(x) = 1 for x G S, and zero otherwise. 
We say that / on R d is piecewise continuous with respect to the decomposition M. d = 
UfcDfc if there exist a finite collection of L disjoint, Lebesgue measurable subsets 
Dk C M d , and bounded, uniformly continuous functions on R d , {fj}f—i, such that 

M d = U^ =1 Dk and / = Ylj=x fj^ Dj ■ If not necessary a part of the terminology such 
as 'with respect . . . ' is omitted. A function fjj on D C M d is piecewise continuous if 
there exists a piecewise continuous / on M d such that fn = f\D. Piecewise constant 
functions are special cases of piecewise continuous ones. 
We consider a 2 nd -order elliptic operator on R d , 

d d d 

A(x) = - d l a t] (x)d J + J2 b 'M d i + O'^ 1 )*) + C ^)' ( J ) 
for which the coefficients must fulfill the following: 

Assumption 2.1 The functions ay = a,-*, b'^ b" (i,j — 1,2, ... ,d) and c are: 

a) piecewise continuous on M. d , c > and a. k j(x) converge to constant values as \x\ 

increases, 

b) there are positive numbers M_, M , < A£ < M, such that the strict ellipticity 

d 

M\x\ 2 < a tJ (x)z l z :j < ~M\x\ 2 , x G M. d (2) 



holds. 



The differential operator Aq(x) = — Ylij=i di a ij( x )dj is called the main part of 
A(x). 

Let us define a real bilinear form on W^(D) x Wp(D), l/p+ 1/q = 1, by: 

d - d - 

a(v,u) = 2~] / aij(x) div(x) djii(x) dx ~2_, / K( x ) v ( x ) diu(x) dx 
, j i • / " i= i Jd 

— I b'l (x) div{x) u(x) dx + / c{x) v(x) u(x) dx. 



(3) 



i=i 



For a domain D with Lipshitz boundary 3D and for each pair v G W^(-D), u £ 
W^(-D) n {Au G (li /1 (Z})) t }, 1 < p < oo, the Green formula must be valid, 

a(u, u) = ( w | A u ). 

The Green formula is also valid for each pair v G W^-D), u G T^ 1 1 (D)n{Au G 7^(D)}. 
The boundary value problem, to be studied in this work, is defined by 

(XI + A(x)) u(x) — fi(x), x e D, . . 

u\dD = 0, { ' 

where A > and D is a domain with Lipshitz boundary. In the case of D = R rf we 
suppose that A > and the boundary condition in (f4]) is omitted. The nonhomoge- 
neous term fj, is a positive Radon measure for p G [1, d/(d — 1)), and /i G (D) for 
P = 2. 

The variational formulation of |4|) for a solution u G (D), p G [1, <i/ (d — 1)) or 
p = 2, has the form: 

\(v \ u) + a(v , u) = for any w G W^(D). (5) 

Solutions of ([4} and ([5]) are called strong and weak solutions, respectively. In the case 
of a problem on K d the variational problem is defined by expression 

X(v\u) + a(v,u) = (v\fi }, for any v G W*(R d ), (6) 

where for p = 2 we have fj, G W^~ 1 (R <1 ) and for p G [1, ^(£>) with a 

bounded D C M d . 

For the differential operator -ff (a;) = XI — a 2 A on R d and A > the fundamental 
solution (x,y) i— > i(A,cc — y) can be represented in terms of the Bessel function 
Ki>, v = [d — 2)/2. Let us denote the corresponding operator by T(A, H), i.e. we have 
(XI — a 2 A)T(X, H) — I on the linear space of continuously differentiable functions 
with compact supports. Much more, for any a > there exists a representation 
of the operator T(A, H) a in terms of an integral operator with a positive kernel 
t a (X, x — y) which is expressed by the Bessel function K v , v = (d — 2a) /2 [Shj . The 
Green function for the differential operator XI — a 2 A on D with the homogeneous 
Dirichlet boundary conditions at dD is denoted by i(A, Hp, •, •) and the corresponding 
integral operator by T(A, Hp). In this way we have (XI — er 2 A)T(A, Hd) = I on the 
linear space of continuously differentiable functions on D. There exist representations 
of T(X,Hu) a as integral operators with kernels t a (Hrj, •, •) which are positive on 
D x D. The equality (XI + Hu)T(X, D) = / enables us to define various closures 
H D = -XI + T(X,H D )- 1 , such as the closure from W%(D) onto W 2 / 1 (D), from 
J)(H) = T(X,H D )L p {D) onto L p (D) etc. 

For A sufficiently large and /i G W / 2 ~ 1 (R d ) solutions to (J3])-© exist and can be 
represented by T(X,H) 1 / 2 (or T(X,Hd) 1 / 2 ) as described in the following. Let us 
define the bounded operator on L2(R d ) by 



W 



T^H) 1 / 2 ^ - M5 l3 )d 3 )T(X,Hfl 2 . (7) 



Then || W||a < (1 — 7) where 7 = M/M < 1, so that there exist the operator 

T(X,A ) = T(X,H)^ 2 (I -W)- l T{\Hf' 2 (8) 
mapping W 2 -1 (K d ) into W}(R d ) with the norm 

\\n\A )\\ L(w -i >wi) < ^(A- 1 +M- 1 ). (9) 

From Aronson's inequalities |Baj we have the following result. The operator T(A, Ao) 
is an integral operator and its kernel is the fundamental solution of differential oper- 
ator XI + Ao(x) on R d . If the lower order differential operators in (QJ are non-trivial, 
then T(X,A) exists for A sufficiently large. Results ((8|) and (J9j) are valid for bounded 
domains as well. We have to replace H with the corresponding Hd and obtain in 
this way the operators Wd, T(X, Ad)- Various closures A, Ad are defined in terms of 
the constructed operators T(X, A), T(X, Aq), respectively, as in the case of differential 
operator — a 2 A. 

In the case of pL G H{D) solutions to ([3])-© also exist. We have the following 
result [BQ1ILR2] : 

THEOREM 2.1 Let D be a bounded domain with Lipshitz boundary. For each p G 
[1, d/{d— 1)) there exists a unique weak solution u of f5[) belonging to the class Wp(D) 
and possessing the following properties: 

(i) There exists a positive number c depending on A£_,M,p,D, such that the following 

inequality is valid: 

\\u\\ P ,i < c/i(D). 

(ii) If {pin : n G N} C 7Z(D) converges weakly to a fx G 1Z(D), then the corresponding 
sequence of weak solutions {u n : n G N} C W^{D), u n — A~^ /i„, converges 
strongly in L p (D) to u = A~f^ pi. 

This theorem is the theoretical background for construction and analysis of con- 
vergence of numerical solutions in Li(JBL d ). 

3 Grids and associated functions 

Let the orthogonal coordinate system in R d be determined by unit vectors , and let 
us define the set G n by: 

d 

G n = {x = h(n) J2 k i e i : k i 5 Z l> ( 10 ) 

1=1 

where h(n) — 2~ n is called the grid-step. A grid-step is usually denoted by h and only 
if necessary by h(n). Elements of G n are called grid-knots and the constructed sets 
G„, n G N are called grids. Sometimes we say that G n discretize M. d . Accordingly, the 
subgrids G n (D) C G n defined by G n {D) = G„ (~l D are called discretizations of D. 
To each v G G n there corresponds a grid-cube C„(l,t>) = J^i [ v j> v j + h)i where Vj 
are coordinates of v G G n . Cubes C„(l, v) define a decomposition of R into disjoint 
sets. Apart from the basic cubes, C n (l,v),v G G n , for constructions we need larger 
sets. Let p G N d . Then 

d 

C n (p, v ) = Y[ [vi,Vi + hpi) 
i=\ 

are apparently rectangles with "lower left" vertices v and edges of size hpi. These 
rectangles define a partition of M d as well. The considered cubes C„(l, v) and rect- 
angles C n (p, v) are semi-closed in the sense that they contain only one of their 2 d 
vertices. 



Basic cubes are defined by their "lower left" vertices. Apart from these basic cubes 
for our constructions we need closed rectangles, 



S n (p,v) = Y[ i v i - h>Pi . v i + h Pi\i (H) 
i=l 

which are defined by central grid-knots v. Apparently, S n (p, v) is the union of closures 
of those basic cubes C n (p,x) which share the grid-knot v. 

The grids G„ of (flQ|) are homogeneous with respect to translations in the direction 
of coordinate axes, i.e x G G n , t = hpiBi => x + t G G n for any i G {1, 2, . . . , d} and 
Pi G Z. There exist subsets of G n which are also homogeneous in the defined sense. 
Let r G G n and r — (n, r 2 , . . . , r^) G N d be fixed. The set 

ri 

G„(r ,r) = i> + h ^ fc, r, e, : k e Z} (12) 
;=i 

is a subset of G„ with the following feature x G G n (ro,r),t = hpiriei =S> x|f £ 
G n (ro,r). A grid (fT2|) is denoted by G n (R), where i? stands shortly for the 2d 
parameters ro,r. 

Let ho = 2~ n ° for some no G N, p G N d and let D be a connected set with the 
structure D = U v& p n C n (p,v), where F n C G n . For the subgrid G n (D) = DCiG n (R) 
the set G n (D) is discrete and therefore its interior, closure and boundary are defined 
indirectly, int(G n (D)) = G n (D)nint(D), cls(G n (D)) = G n (R)nD and bnd(G n (D)) 
is the difference of cls(G n (Djj and irvt[G n {D )). Apparently, int(G n {D)} C G n {D) Q 
cls(G n (D)) . Let a finite collection of sets Z?;,Z G £ make a partition of K d , where 
each Di has the structure like the described set D. Then G(l) — Di n G n {R) make a 
partition of G„. 

Each a; G G n can be indexed by m G M d , where a? = hm. Similarly, we index 
grid-knots of G„(ro, r) by those m G Z d for which there holds cc = ro + m i r i e i- 
Therefore, we define the sets I n = TL d and I n {R) G /„, indexing the grid-knots of G n 
and G n (R). In this work frequently utilized pairs of grids and their index sets are 
Gn , In , G n (R), I n (R); G„{1), /„(/); G n (R,D),I n (R,D). 

The s/ii/t operator Z(x),x G M d , acting on functions / : R d h- » K, is defined 
by i y Z{x)f s )(x) = f(x + z). Similarly we define the discretized shift operator by 
(Z n (r, i)u n ) k = (u„)j, where I = k + rhei. 

Discretization of differential operators. With respect to a grid step h, the 
partial derivatives of u G C^($L d ) are discretized by forward/backward finite differ- 
ence operators in the usual way, 

di(t)u(x) = \(u{x + tei) - u(x)), d . . 

a«(t)«(x) = \{u{x) - u(x - tet)), X £ K ' f t U ' {l6) 

Letr G Z\{0}. Discretizations of the functions diU on G„, denoted by [/j(r)u„, T^(r)u n , 
are defined by: 

(Ui(r)u n ) m = Bi(rh)u(x m ), (V l {r)u n ) m = Bi(rh) u(x m ), 

where x G G n . Then 

Ui(r) = (rh)-\Z n (r,i) - I), , . 

Vi(r) = {rh)- l (l-Z n {-r,i)) = U t (-r) = -U t (r) T . 

Therefore we have Ui(—r) = f/j(r) Z n (—r, i) = Z n (—r, i) Ui(r), and similarly for Vi(r). 
In the case of r = 1 we use a short notations E/j, V^. 

In accordance with the previous terminology, we say that <9j, ■ dia^dj etc. are 

differential operators on R d or Z?. We say that their discretizations are defined on G n 
or G n {D). In particular, discretizations of the differential operator ([T]) are denoted 
by A n . Naturally, matrices A n are the main object in this work. 



3.1 Relations between l(G n (R)) and W^-spaces 



The discretization of a function u G C(M d ) on G n is denoted by u„ and defined by 
values at grid-knots, (u„) m = u(x m ) where x m = (mi/i,m 2 /i, . . . , m^/i) G G„, and 
m = (mi,m2, . . . , md) is a multi- index. The function u„ is usually called a grid- 
function. We denote the linear spaces of grid-functions by l{G n ) or l(G n {D)). The 
linear space of grid- functions on G n (R) with finite supports is denoted by lo{G n (R)). 
Elements of l{G n ) are also called columns. The corresponding L p -spaces are denoted 
by l p (G n ) or l p (G n (D)), and their norms by | • \ p . The duality pairing of v G l q {G n ) 
and u G l p {G n ) is denoted by ( v|u). The scalar product in h{G n ) is denoted by ( - |- } 
and sometimes by (-|-). The norm of l p (G n (R)) is denoted by | \r p . For p G [1, oo) 
this norm is defined by: 



\u\ Rp = vol(R) £ 
kei n (R) 



i/p 



where vol(f?) = Y\i=i r i- Finally, for p = oo we have |u|ftoo = sup{|ufe| : k G I n (R)}. 

Let us define the quadratic functional on l(G n ) by q(u) — | u 1 1 and 9k(u) = 
vo\{R)Y?i \Ui(ri)u\ 2 R2 on l(G n (R)). It is understood q R = q for G n (R) = G n . A 
discrete analog of W^-spaces is the linear spaces of those u„ G l(G n (R)) for which 
the norm | • \r,2,i- 

|u|ki = I u lfi2 + «fl(u), (15) 

is finite. This space is denoted by w\{G n (R)). By convention | • | 2 ,i = | • \r2,i for r, = 
1. The subspace of grid-functions u G w\(G n (R)) for which u„ = l Gri ( D }\i n is denoted 
by wl(G n (R,D)). Hence, w\{G n (R,D)) for r = 1 is denoted by «4(G„ (£>)). For 
problems on bounded domains we need a discrete version of the Poincare inequality 
which is formulated as follows: 

LEMMA 3.1 Let D be bounded. Then the norms \ ■ 1 2.1 and gij(-) 1 ^ 2 are equivalent 
mw\(G n {R,D)), 

te(-) 1/2 > P\-\r2,i, 

where (3 is independent of n. 

An element (column) u„ G l(G n ) can be associated to a continuous function on 
M. d in various ways. Here is utilized a mapping l(G n ) 1— ► C(M. d ) which is defined in 
terms of hat functions. Let \ be the canonical hat function on K, centered at the 
origin and having the support [—1,1]. Then z ^ (f>(h,x,z) = \{h^ 1 {z — hx)) is the 
hat function on M, centered at x G R with support [x — h, x + h] . The functions 
z l— > 4>k{z) = rii=i 4>{hi x ii z i), x i = ^^ij define c?-dimensional hat functions with 
supports 5 n (l, a;) = Jlil^i — h,Xi + h}. The functions 0fc(-) G G„, span a linear space, 
denoted by E n (M. d ). Let u„ G l{G n ) have the entries u n k = (u„)fe. Then the function 
u(n) — Yl,k&i n u -nk(t>k belongs to E n (M. d ) and defines imbedding of grid-functions 
into the space of continuous functions. We denote the corresponding mapping by 
$„ : l(G n ) 1 — ► E n (R d ). Obviously that there exists : E n (M. d ) ^ l(G n ) and 

the spaces l(G n ) and E n (M. ) are isomorphic with respect to the pair of mappings 
$„, It is clear that E n (R d ) C £„ + i(R d ) and the space of functions U n E n (R d ) is 
dense in L p (R d ),p G [1, 00), as well as in G(R ). Let us mention that J2k '/'fc = 1 on 
R d . 

Now we consider another collection of basis functions. To each x = hk G G n (R) 
there is associated a <i-dimcnsional hat function 



= n 



X 



- hkj 
hr; 



obviously, with the support S n (r, x) =\\ i [ + rih]. They span a linear space 

denoted by E n (R,R d ). Again we have J2k^k = 1 on R d . The mappings ^n,^ 1 



cannot be applied to elements of l(G n (R)) and E n (R,M. d ), respectively. Therefore 
we define restrictions & n (R) : l(G n (R)) — ► E n (R,R d ) and the corresponding inverse 
mapping $~ (i?) by the following expression: 



u(n) = $„(i?)u„ = ^2 ( U n) k ^k- 



(16) 



If we have to underline that u(n) is related to a particular set of parameters R 
then we use an extended denotation u{R,n). For two functions v(n),u(n) we have 
(v(n)\u(n)) = h d vol(R)J2ki s ki v kUk where s k i = \\ipk\li 1 (ipkl^i)- Lct us notice that 
Ez s ki = !• 

We cannot compare directly columns u„ with various n. An indirect comparison 
can be made by using u(n) = <I> n (i?)u n £ U n E n (R, R d ). To compare Ui(ri)u n and 
diu(n) we need an additional expression. Let u and u(n) be related by (|16p and 
Ski(i) — h" d (di^pj e \di^i). Then, for a homogeneous grid G n (R),r e N d , there must 
hold 



V kUlS k l(i) = -| ^ 



Uk+r') S 0r i So r 



(17) 



where r' = (n, r 2 , . . . , r-j-i, 0, r i+ i, . . . , r d ), s 0r = (d i i>o\d i ip r ) and s' Qr , is the (d - 
l)-dimensional quantity sm- Then (| 1 T[) follows from Ez = 1 an d consequently 
Ez *fcz = 0j after the sum is carried out over any partial component li of the index I. 
Thus we have 

{d l v{n)\d l u(n)) = ||Vfe]U E s fe'*' ( U ^»m ik ' i U ^) u ) mi l" ( 18 ) 

where the indices are defined by k' = (fci, . . . , fcj-i, fei+i, • ■ ■ , fed) and analogously i'. 
LEMMA 3.2 Let sequences of functions v(n),u(n),n 6 N 7 be defined by 116)] . XTien 



(0 



< h d 



(ftv(n)l^(n)) - ^vol(iZ) E fee / n ( fl ) (^(n)v) fe {U j (r j )u) k 

\U l {r i )v\ R2 sup f\(Z(w,j)-I)Uj(rj)u\m ■ \w\ < rjh} 
\Uj(rj)n\ia sup { \(Z(w, i) - 1) Ui{n)v\ R 2 : \w\<r 2 h}. 



Proof: After applying the CSB-inequality to l|T8j) and using J2i s ki = 1 we get 
(i). Assertion (ii) is proved here for i = j = 1. In this proof d = d\. By (fT8|) we can 
straightforwardly calculate 



| du(n)) = ||V>fe||i ^ s 0r > (Uj(ri)\) h (Ujjr^u) 



k.r' 



where 



5(n) 



ll^lli £ sow (Ui(r<)v) fc [(^(rOu)^ - (^(n)^ 



By the CBS inequality the error term 5{n) can be estimated as in Assertion (ii).QED 

LEMMA 3.3 Let G n {R) be a homogeneous subgrid defined by ilty) . There exists 
a 2 e (0, 1) such that 

(1 - a 2 )h d \u n \ R2 < \\u(n)\\l < h d \u n \l 2 

uniformly with respect to n 6 N. 



Proof: Let us consider first the one-dimensional case. The grid G n (R) consists 
of points Xk = hrk £ 1,1; e Z, and E n (R,M) is spanned by the hat functions ipk 
centered at xj. with the supports [—hr + + hr] . We define the matrix 5(1) with 

entries: 

1 f (2/3) for fc = /, 

«« = Vr {M^) = ( (1/6) for k = l±r _ 

Obviously we have 5(1) = I — (1/3) A, where the matrix A has the structure A = 
1+ (1/2)(J_|_ + J_) and I± are the first upper and lower off-diagonals. It is well known 
[Str that A has a purely continuous spectrum in [0, 2] so that 5(1) has the spectrum 
equal [1/3,1]. Therefore 



\u{n)\\\ = hvol(R) ^2 Skiu k ui > ^hvol(R) ^ ul = ^h\u\ 



3 

k 



2 

n-2- 



Hence, we have here 1 — a 2 = 1/3. 

In order to generalize this proof to d-dimensional case we proceed as follows. The 
symmetric matrix S(d) with entries Ski can be represented as the outer product of 
d matrices 5(1) with entries as in the first part of proof. Therefore its spectrum is 
Sp(S(d)) = n,-=i Sp(SW)- According to the first part of proof the matrix 5(1) has 
its spectrum in the interval [1/3, 1], implying Sp(d) > 3~ d . Hence, with a 2 = 1 — 3~ d 
we have 

Vol(i?) J2 SklUkUl > (1 - cr 2 ) |u|^ 2 , (19) 

providing us with a proof of the general case. QED 

THEOREM 3.1 Let u(n) = $ n (R)u n . There exists a 2 e (0,1), independent of n, 
such that 

(l-a 2 )h d \u n \ 2 R2>1 < MR^Wl, < h d \u n \ 2 R2)V 

Proof: It is sufficient to prove the first double inequality. The estimates from 
above are obvious. To get the estimates from bellow it suffices to consider diU. From 
Expression (fT5|) we have 

\\diu{n)\\l = U k h J2 E s *'*' ( U i( r i) u )mv ( U ^) u ) m i'- 

m k'V 

Then after applying (I19p to the inner sum we get 

||S W (n)||| > (l-OH^Hr £ {UM^L, = (l-a^h^U^ull,, 

m,h' 

from where follows the estimate from bellow. QED 
An element u G W 2 1 (M d ) does not necessary belong to E n (R, M. d ). In order to 
approximate u with elements of E n (R,R d ) we define: 

u{n) = £ {^k\u)^ k - (20) 

kei n (R) 

The numbers ||'!/'fc|lr 1 (V'fcl u ) are called the Fourier coefficients of u. 

The basic result for our proof of convergence of approximate solutions is formulated 
in terms of the quantity T p (w, u) defined by: 

r p (ii>,u) = \\(z(w) - i) u \\ p 

The kernels 



define an integral operator which is denoted by K n . Actually, the kernels io n define a 
5-sequence of functions on M. d x M. d and K n converge strongly in L p -spaces to unity: 



COROLLARY 3.1 Letpe[l,oo\. Then 

(i) \\K n \[ p < 1. 

(ii) There is a positive number n(R) 7 independent of n, such that — K n )u\\ p < 
k(R) swp{T p (w, u) : \wi\ < hr{\. 

(Hi) The operator K n £ L(L 2 (M. d ), L 2 (M. d )) has the spectrum equal Sp{K n ) = {0} U 
[3- d ,l]. 

Proof: Only (iii) has to be proved. The symmetric operator K n is reduced by 
E n (R, R d ) and represented by an integral operator with the kernel (f2"T1) . It is zero oper- 
ator in the orthogonal complement E n (R, M d )- L . With respect to the mapping < &„(i?) 
the operator K n is mapped to the symmetric matrix K n = ^ n {R)~ 1 K n ^ n (R) = S(d) 
mh{I n {R)). QED 

THEOREM 3.2 Let v,u £ W%(R d ) andu(n),v(n) be defined by Then 
\(v(n M n)) - (v\u)\ < c(R) min { | ^ Z^^vtu)' 

{ || V \\2 SUP\ w \< hr \ 1 2{W, U), 

\(div(n)\dju(n)) - (div\dju)\ 

\\div\\ 2 \\djU - 3ju\\ 2 + sup T 2 (w,dju) 

^ iT}\ ■ ) L \w\<hr\ 

< c(R) mm ^ r i i- i 

\\dju\\ 2 \\diV - BiV\\ 2 + sup r 2 (w,div) 



\w\<hr\ 



where c{R) is n-independent. 



A proof of this theorem is rather technical. For instance, in order to prove the 
second inequality one has to use a sequence of replacements: (div(n)\dju(n)) — > 
(3i(h)v\K n 3j(h)u) — » (3i(h)v\K n dju) — ► (div\K n dju) — > (div\dju). Each replace- 
ment gives rise to an error. The sum of errors can be estimated by an expression as 
given in the second inequality of assertion. 

3.2 Imbedding of l(G n (R)) into W / 2 ~ 1 -spaces 

Beside the functions i/'fc we consider the functions defined by: 

Xki+{x) = i[kiM+nh){Xt) Y[ ^{Xj), Xki-{x) = Xk-n<Hi+( x ) 

for all the possible i = 1,2, ... ,d, and the linear space F n (R,M. d ) spanned by the 
defined functions Xki±- The integral operators Kn\xiX) with the respective kernels 

u$(x,y) = V] n — — n- Xk l+ {x) xkt+{y) = V I, — ^ ^ Xkt-(x)xM-(y) 

have properties similar to the integral operators K n of Corollarv l3.ll The same is valid 
for the non-symmetric integral operators Kn\ip,x)i an( i their adjoints Kn\x,ip), 
where the kernel ofKn(ip,x) is defined by (x,y) i-> (ft d vol(-R)) _1 J2k i ) k{x)Xki+{y) ■ 

LEMMA 3.4 The operators K$ (x, x), Kn^ (ip, K n l \tp, xV have properties (i) and 
(ii) of Corollarv \3.1[ The spectra of operators K n (XjX) S \-(L 2 (M. d ), L 2 (M. d )) consist 
of and an interval [re, 1] with certain re £ (0, 1). 

An element fj, £ W 2 ~ 1 (M d ) is represented as /x = f + d ifi witn fi e L 2(^ d ) 

with the norm |[/i[[2,-i = (X^ ll/illl) 1 ^ 2 - Its discretizations are defined by grid- 
functions fi n with the components: 

« = B ^(««- R bt <a ' w,) - (22) 



Therefore we can write: 



(u(n)\(i) = h d (u„ | fi n ) R = \\ip k \\i ^2 MfeMfe- (23) 

feeJ„(R) 



The expression 

diipk = — r 

Tift 



Xki- — Xki+ 



a.e. on K d , 



enables us to rewrite the components of fj, n in (|22|) in terms of the Fourier coefficients 
fik± = HV'fellr (Xki±\fi) 01 functions fi. The corresponding grid-functions f in ± arc 
imbedded into the space F n (R, K d ) by the mappings fj n ± i— > fi±(n) — Kn\x,X)fi 
which are completely analogous to the mapping fo fo(n) — K n fo of Corollary 13. II 
Now we can get the following useful expression: 

d 

(u(ri)\n) = h d (u„|f 0n )H - ^2 (Ui(r i )u n \f in+ ) . (24) 

2=1 

Therefore, by (|14[) we can represent elements fi G l(G n (R)) in the following way: 

d 

= fori - ^ ^( r i) f i"+- ( 25 ) 

i=l 

From Expressions (23) and (M) we have |(u„|/n n )| < |u n | 2 ,i(|fo|| + |W|1) 1/2 , 

implying a natural definition: 



Kll-i = El f "l2 + lwl 



For the sake of a concise writing of final results we denote here K n by Kn^ 1 ■ 

LEMMA 3.5 For each fi G W^QR ) and the corresponding discretizations \i n de- 
fined by \25)) the following is valid: 

and 

o < ||mII1-i - ^iMnlU = J2(fi\(i-K%Hx,xM) - o, 

i=0 

as n — > oo. 

3.3 Imbedding into Wp-spaces 

Theorem 13 . 1 1 can be partially generalized. 

LEMMA 3.6 The following assertions are valid: 

(i) If p G [1, co] and u„ G l p (G n (R)) then: 

\\u(n)\\ pA < h d ^\u n \ RpA . 

(ii) If p G [1, co] and u G W^(R d ) then: 

\\u{n)\\ p ,i < h d ^\u n \ RpA < \\u\\ pa . 



(in) Let u = u + +u , where u + — max{u,0},« = min{u, 0}. Then 

h d \n n \ R1 < || u ||i = Wu+i^Wx + Wu-in)^ < h d (\u+\ m + \K\rn)- 

Proof: Let us consider the mapping : l{G n {R)) i-> E n (R,R d ) denned by 

u(n) = Ukipk- In terms of functions Xki+ of Subsection 13 . 21 we easily get: 

d l u{n) = ( U i(Pi) u n) k Xki+- (26) 

kel n (R) 

Apparently we have for p = 1 and p = oo the following inequalities: ||u|| p < 
/i d / p |u n |fl p , ||<9ju|| p < h d / p \Ui(pi)VL n \ Rp . Hence, by the Riesz-Thorrin theorem we 
get (i). 

Similarly we prove (ii). For p = 1 and p = oo we have /i d/,p | u n | Hp < II u Hp- Now 
we use the expression ipk +pei (x) = ipk{x ~ phei) and get 

(Ui(p)u n ) = (ip k | Bi(p)u). 

Wk\\l 

Because of \\3 l (p)u\\ p < \\diu\\ p for p = l,oo, we have h d / p \Ui(p)u n \ Rp < \\diu\\ p . 
Assertion (ii) follows now from the Riesz-Thorrin theorem. 

The right hand side of double inequality (iii) is implied by (i) while the left hand 
side is implied by (ii). Therefore we have to prove the equality ||w||i = ||u + (n)||i + 
|u~(n)||i. This equality is a consequence of J2k^ k = 1 on R d and = 
E k (^k\u+). " QED 



4 Construction of discretizations 

It is important to underline at the beginning that discretizations A n of differential 
operator A(x) are defined prior to discretizations of the forms ([5]), ([6]). This fact 
is in a full agreement with construction of discretizations A n in this section. Some 
classes of discretizations are derived from a general principle which is not based on 
finite difference formulas and cannot be apriori related to variational equalities. Nev- 
ertheless, bilinear forms must be associated to A n so that A n are derived from the 
corresponding variational equalities. The constructed bilinear forms are considered 
as discretizations of the original form ([3]). These forms are basic objects in our proof 
of convergence of approximate solutions. In the next two subsections schemes and 
the corresponding discretized forms are constructed for two classes of methods. 

Forms a„(-, ■) on E(R,WL ) x E(R,R d ) and matrices A n on G n (R) are related by 
equalities: 

a n (v,u) = (v\A n u) R . 
In addition, the discretized forms determine the discretized variational equalities: 

A(v„|u„)_r + a n (v,u) = (v n \n n ) R , v„ S w 2 ,i(G n (R)). (27) 

To discretize A(x) means to associate to A(x) a sequence of matrices A n on 
G n (R),n £ N. Of course, the matrices A n must be constructed reasonably in order 
to enable demonstrations of the convergence of numerical solutions. The convergence 
analysis is postponed until two next sections. Therefore, in this section, the termi- 
nology "discretizations" of A(x) instead of approximations of A{x) seams to be more 
suitable. 

Discretizations to be considered in this section are possible if certain conditions 
on a,ij are fulfilled. The required conditions are stronger than in Assumption ^. II By 
relaxing them gradually as n — > oo we obtain discretizations for a general A(x) given 
by Assumption 12.11 

To a given diffusion tensor a = {aij}f d we associate an auxiliary tensor a defined 
by the expressions: 

an = an, = —\aij\ i^j. (28) 



Assumption 4.1 (on discretization conditions) 

1. There exist q G N d , a finite index set £j and a partition M. d = U;-D;, where the 
sets Di are connected unions of C n (q,x),x E G n (Q). The sets Di and the diffusion 
tensor a = {cii/} ff must be mutually related as follows: 

- The tensor-valued function x t— > a(x) is continuous on Di and the functions 
a ij ? * 7^ 3 do not change signs on Di . 

- There is a grid-step h e = 2~™( £ ) and the closed sets 

Di(e) = UxedsiD,) S n ( £ )(q,x), 

such that the functions ciij can be extended to Di{e), not changing the signs 
on Di(e), and the strict ellipticity is valid on Di(e) with the same bounds 
M,M. 

2. For each I G £ the auxiliary diffusion tensor a is strictly positive definite on Di (e) . 

3. To each Di there is associated parameter p(l) G N d . such that the following in- 
equality is valid: 

uj(a) = inf n minj * inf { — — inf au(z) 

I Pi(l) *£S n (p(l),m) 

- ^ — ttt sup \a im (z)\ : x e G n n Di(e)\ > 0. 

Pm{l) zeS n {p{l),x.) J 

Condition 3. is crucial in our construction of discretizations A n which have a 
particular feature called the compartmental structure. In the next definition / is an 
index set, and matrices A = {a^}// are considered in linear spaces l p {I) consisting of 
functions on the set /: 

DEFINITION 4.1 (Compartmental structure) A matrix A in ioo(-0 « s said to 
be of positive type if A = pi — B, p > 0, B > and |-B|oo < P- It is called 
conservative if Bl = pi. A matrix A in h(I) is said to have the compartmental 
structure if A = pi — B , B > and \B\i < p. It is called conservative if for each 
u > there holds |-Bu|i = p|u|i. 

LEMMA 4.1 Let A in l\(I) be compartmental. Then 

(i) If A is conservative then |-B m |i — p rn for each m G N. 

(ii) The spectrum sp(A) of a compartmental matrix is contained in 5fiA > 0. If A is 
conservative, then G sp{A). 

In the case of d = 2 Conditions 2. and 3. can be always fulfilled. For d > 2 there 
exist positive definite diffusion tensors {a,j}ff such that Conditions 2. and 3. are not 
possible |LR3] . 

Discretizations A n are defined in terms of its matrix entries (A n )ki, where hk, hi G 
G n . For a fixed x = hk G G n the set of all the grid-knots y — hi such that (A n )ki ^ 
is denoted by J\f(x) and called the numerical neighbourhood of A n at x G G n . In our 
constructions the sets J\f(x) for x G G n PiDi are mutually alike. A set N(x) contains 
always a "cross" consisting of x and 2d elements ±hpi(l)ei. Additional elements of 
J\f(x) depend on the sign of ay , i ^ j. 

Two classes of discretizations are analyzed. One of these classes can be defined 
straightforwardly in terms of forward and backward difference formulas. The resulting 
discretizations are called basic schemes. Discretizations of the other class come from 
a general principle [LR3] and they are called extended schemes. 

The compartmental structure of discretizations A n of the differential operator 
A(x) is the goal of overall analysis. Here we describe a general approach to the 
constructions of discretizations A n with the compartmental structure which is based 
on reduction to the two-dimensional problems. 



Figure 1: Asumption on function 012 



oia > 



a 12 < 



oia < 



oia = on UC n (p,v) 
Pi = 2, P2 = 1 



The index set of pairs 1(d) = {{ij} : i < j, i,j = 1,2, ... ,d,i ^ j} has the cardinal 
number m(d) = d(d — l)/2. To each index {/ci} G 7(d) we associate three coefficients, 

a *;fc Y a fcfc, a ;i ~~ Y a "' a M _afci > ( 2y ) 

and a bilinear form o^ fc ^(-, •), 

a {kl} [v,u) = ^ / ajf } (aO<9M*)^(aOda:. (30) 

*,jS{r,a} D 

Apparently, for each pair u, u S C^QR. ) with compact supports, the following equal- 
ity is valid: 

a(v,u) = a {fei} («,«). 

{fci}ei(d) 

To each of the forms a^ kl ^(-, •) we must associate a sequence of forms On^(-, •) and 
matrices -An^ constructed by two-dimensional schemes. Then the matrix 

A n = 4 W} . (31) 

{fei}e/ 

is a discretization of j4o(a;). If each An kl ^ has the compartmental structure then A n 
also has the compartmental structure. However, A n can have the compartmental 
structure although no An kl ^ is compartmental. Condition 3. of Assumption 14.11 
ensures this advantageous property in our constructions. 



4.1 Two methods of discretizations 

The forms must be constructed by the rules (|29 p -(|30 p having in mind that the con- 
struction for higher dimensional cases depends on the construction for two-dimensional 
case. Therefore, we are due to describe the construction for the two-dimensional case. 



Basic schemes 

In the case of numerical grids G n (P) the numerical schemes and corresponding dis- 
crete bilinear forms can be easily mutually related. A proof of convergence in Sobolev 
spaces for various right hand sides are presented in our works [LR2, LR3]. In the 
present case we extend analysis to problems with the function 012 having both signs. 
Difficulties appear at those grid-knots where the sign of a\2 changes. 

The discretizations are constructed by assuming that the sets {x G R d : 012(0;) < 
0} and {x £ R d : 012(2;) > 0} are separated by a set (connected) so that it contains a 



connected subset which is equal to the union of cubes C n (p,x),x E G n (P) for some 
n G N. An illustration of this assumption is given in Figure [1] This assumption 
is not valid for a general matrix- valued function x i— > a(x). Therefore we have to 
comprehend this assumption as a step of an approximation procedure in our process 
of construction of discretizations A n . For the sake of simple and brief presentation 
we assume in the next construction that the assumption is valid already for n = 1. 
Let us define the sets D n (—) C K d by: 

Ai(-) = {LWg„(p)C„(p,'u) : ai 2 (x) < for x e C n (p,v)}, 

and D n (+) — D n (—) c . Now we define the subgrids G n (P,—) consisting of all the 
vertices v E G n (P) which determine the set £)„(—). As well we need the subgrid 
G n (P,+) consisting of those grid-knots x for which the segments I(x,y),y = x + 
hnPi^i have the following property I(x,y) C cls(D n (+)). It is easy to verify that 
each x G G n (P) must be contained in one of sets G n (P, =p) and each of segments 
Ii(x,y),y = x + hpiGi must have both of its end points x,y in some G n (P t =f). 
Some of grid- knots and some of segments are contained in both sets, G n (P, =f) and 
cis(-D n (=F)), respectively. For a segment in e2-direction, I(xi, X2), x% — X\ + h n p2e2, 
the following is true. If X\ G G n (P, — ) then /(cci, #2) C cls{D n {—)). If a; 2 G G n (P, +) 
then 7(^2, cci) (downward vertical segment) may be outside of cls(D n (+))). If this 
happens then this segment is contained in a cube cls(C n (p, z)) on which the function 
(Z12 has zero values. This fact is a consequence of a strict separation of supports of 
functions max{ai2,0} and min{ai2,0} and will be utilized in our constructions of 
discrete bilinear forms. 

To define forms and entries we need 

z (±) = Piei ± pjej G /„, 

x (±+)( n ) = l(±hp l e l + hpjBj) GS n (p,0), 

a;( ±_ )(n) = |(±/ipiej - hpjej) £ S n (p,0). 

Obviously we have hz^ a ^ G G n while x^\n) are not necessary in G n . 

Let us consider a sequence of two-dimensional forms on E n (R,R 2 ) x E n (R, R ) 
which are defined by the following expressions: 

a n (v,u) = ah\v,u) + a,n(v,u) 



a 



a 



{v,u) = Eij=i Exgg„(p,-) {3i(pih)v)(x) 

x a tJ {x + x^ ++ Hji))(B 3 { P] h)u){x), (32) 

(+)(«,„) = £*. =1 E» eG „(p,+) (a,((-irV,%)W 

Discretizations A n of differential operator Ao(a;) can be easily obtained from the 
constructed forms variationaly. Let us define matrices A„ (i,j,Pi,Pj) as the diagonal 
matrices with entries A„ (i, j,Pi,Pj) X x = a<ij{x + x^ + ^ (n)). In terms of matrices 
Ui(pi),Vi(pi) and Z n (pi,i) of (fT4|) we get the following expressions: 

2 

= - J2 J2 ViiPJ^nHhj'PhP^UjiPj)^^,-) 

leL- y=i (33) 

- 51 Z F ^) A ^ A{ nHi,3,Pi,Pj)A-jU j (p j )l Gn (P,+), 

where Ai = I, A 2 = Z n (— p 2 >2) and 1g„(p.=f) are the projectors on the linear sub- 
spaces of grid-functions with supports in the sets G n (P, =p), respectively. The entries 
of A„ can be also easily calculated for grid-knots of int(G n (P, =f)). In order to get 
simple expressions we use the following abbreviations: 

a ( f ] = aij (x + x^{n)), a,0e {+,-}. 



Figure 2: Numerical neighbourhoods at two grid-knots of G n (P, =f) 




012 > 




Pi = 3, P2 = 1 



012 < 



The entries on the cross x ± h r p r e r , r — 1,2 have the structure: 



(A 71 

(A r , 



kk±pie± 



kk±p2G2 




for <zi2 < 0, 
for ai2 > 0, 



The entries in the plane spanned by e%, e2 have the structure: 



(Avfefe=bs( 

(A- 



TT~ — krf^l f° r °i2 < on Di, 
h 2 pip 2 1 1 



fcfe±z<+) 



h 2 PiP2 



l 12 



for ai2 >0onfli, 



where is associated with the upper and —z^> with the lower indices of a^ a ^ . 

The diagonal entries (A n )kk are equal to the negative sum of all the entries (A n )ki, 
I k. Entries for grid-knots at bnd(G n {P,T)) can be more complex. In two- 
dimensional problems the forms (|32p are natural and the entries of A n calculated from 
variational equalities (f2"T)) can be perceived as buildups made by forward/backward 
finite difference formulas. This approach is self-understanding and we call it the stan- 
dard approach. The result are basic schemes of discretization. An illustration of 
numerical neighbourhoods of basic schemes is given in Figure These schemes are 
studied thoroughly in SMMM . The strict ellipticity of forms a„(-, •) is expressed 
always in terms of the original pair of positive numbers M_, M . Generally, the strict 
ellipticity of discretized forms follow from the compartmental structure of A n . These 
problems are analyzed in the next section. 



Extended schemes 

Contrary to the standard approach in developing discretizations of Aq(x) we have 
methods based on some general principle and which are not a priori related to the 
forward/backward finite difference formulas. One of such methods is described here. 
Principles of construction are given in [LR3] . To each Di we must associate elements 
r(l) £ N playing the role analogous to p for basic schemes. 

Again we assume the strict separation of sets F(±) = supp (max{±ai2, 0}). The 
set L of Assumption 14. II is partitioned into the subsets L T , where I G £_ means that 
Ojj < on Di and I G £+ means <Zjj > on Dj. Let us remind that the sets Di in 
present case are determined in terms of cubes C n t e -\(r{l), x). There is always a room 
of arbitrariness in a determination of these sets. The following maximal property 
removes some of arbitrariness. There exist a n(e) G N such that the sets Di,l G £_ 
have the following properties: 

a) ai2 < on Ui<=c_Di, 



b) Each Di satisfies Assumption 14.11 



c) Each C n r e )(r(l),x),l E £_,cc € G n on which 012 = is contained in some of 

Di,leL- 

Now we can define D{—) = Uien_D[ and D(+) = D(—) c . In our proceeding dis- 
cussion we assume n(e) = 1. The subgrids G„(— ) = G n n D(-) have the same 
properties as the corresponding subrids G n (P, —),pi = P2 = 1 in the subsection on 
basic schemes. Similarly, the sets G n (+) = G n nD(+) coincide with G n (P, +) as well. 
Now we define subgrids G n (l, — ) = G n (— )C\Di and conclude that G n (l, — ) form a par- 
tition of G n (— ). However, the sets G n (l, +) = G„(+) fl Di do not form a partition of 
G n (+) because some of grid-knots at dD(+) may be outside of each G n (+, I). There- 
fore, we have to extend the sets Di to wider sets Di such that G„(l, +) = G n (+) DDi 
form a partition of G n (+). The sets Di cannot be defined uniquely. Here we demand 
the following properties. The sets must be disjoint, and the closure of int{D{) must 
coincide with cls(Di). Thus we have: 

G n (i,-) = G„(-)nA, G n (i,+) = G n (+)nA. 

In accordance with our discussion about properties of sets G n (P, ±) in the subsection 
on basic schemes we finally conclude that G n (l, ±) cover G n and some of them may 
have common grid-knots. 

The forms «„(•,•) are expressed in terms of a,ij(x + x( '(l,ri)),x £ G„, where 
x( ±:iz \l,n) are certain elements in R d . Since x^ ±±s> {l,ri) can take values outside of 
cls(Di) we are due to specify how to take values of aij [x + a;( ±± ' (I, n)) in such cases. 
The values must be taken in the set Di(e) of Assumption ^. II In this way we conclude 
that the entries of A n are calculated in terms of values of coefficients at points 
which are not grid-knots. 

It is convenient to use a representation a n {v,u) — an, \v,u) + a„ where 
the forms a^\v, u) are related to the index sets as previously. Let us define 

t(±+)(r) = \{±hr l {l)e l + hr {l)e ) G 5„(r,0), 
t(±-)(r) = f(±hr i (l)e i - h rj (!)ej) € S n (r,0). 

Obviously x^ a ^ (I, n) and t^ al3 \r) coincide for p = r(l). The form a„ '(■, •) is defined 
by: 

2 

ie£_ xeG„(0 »=i 

+ + t {++) (r)) (a t (nm)v)(x)(a 3 (r 3 m)u)(x) 

i*i (34) 

+ ]T ^ (a . + t(++)( r ))IiM [(a i (%)(x)(a,(/ l ) u )(x) 



- (a i (r i (/)/i)«)(x)(a i (r i (0%)(x) 

We obtain al + ^(w,u) from an (v, u) by replacing Bi(ri(l)h) with 3j((— iy~ 1 (r i (l)h)), 
£_ with £ + and t( ++ ' with t^" 1 The forms a^\v,u) are not second degree poly- 
nomials of 3i(h) with simple structure. Due to the compartmental structure of A n 
they can be ultimately represented as forms depending on Bi(qih) with various 
For the quantities a^(u,ii) more comprehensible expressions can be written down 



such as dHjl and (|48j) . 

Discretizations A n have a general expression: 



where 



i=l i^j 

£ - ^(r i (I))^+)(i,i,r i) r J -)Di(r < (I)) 



A { rt\l) = - ELi ^A i (r i ) r ^- ) (t,i,l,l)A i (r i )^ 



J n(l)h 



- £ ^(r i (0)A i (r i ) T 4-)(i,i,r i ,r j )A j (r j )f/ J -(r j (/))- E , „„, 

x [^(^^(^M)^)^ - V i (r i (0)^- ) (i.i,»- i ,r i )Ui(r i (0) 

and where Ai(ri) = /, A 2 (r 2 ) = Z n (— r 2 (Z), 2) as in the previous case. 

In order to write down the entries of A n we need the following abbreviations: 

- rj(l)ej £ I n , 
t< a «(r)), a,p G {+,-}, 
-t(++)(r) - /iei), 
- t (++)(r)-/ie 2 ), 
-), a|-- } (r) = a£r-) (r .). 

Then we have the following nontrivial off-diagonal entries of A n : 



(35) 





= n(i) 




— O-ij ( 


«12 + V) 


= ai2 




= a\2 






= a it 



{A r , 



^ +) (D-?il^ +) (r)| for « 12 <0, 



i 



r 2 (/) 

,(+±) m r a (0 u(+±) 



1 " I 1 ) 77T a i2 ( r ) for ai 2 > 0, 36 

r 2 (i) 



«U ±e2 = -^[4^(1) - ^j|air ; (r)| 

The entries in the plane spanned by e\, e 2 have the structure: 

{ A n) kk±w( - ){1) = ~ Jfrjtfcffi \ a i2 T) ( r )\ for ai2<0onA, 



(37) 



where +u/ ± ) is associated with the upper and -ujW with the lower indices of a'"^, 
respectively. 

Discretizations A„ defined by (|36|) . (|37| are called extended schemes. The numer- 
ical neighbourhoods are illustrated in Figure [31 

In (|35[) we have 4 sums with respect to the indices The first and second 

sums have expressions similar to Expressions (f3"2"j). They contribute to Expressions 
(|37[) and to a part of entries in (|36p . Unfortunately there appear non-trivial entries 
(^4n)fci, £ = fc ± Ti&i. These entries must be canceled by contributions from the third 
and fourth sums. These two latter sums lack the structure similar to (I32p since the sum 
includes the terms a^BivBiit and ai 2 3 2 w3 2 w. So the net result of all four sums are 
entries ([37]) . Let J C I n {l) consists of grid-knots in some of x £ G n (l, ±) for which the 
numerical neighbourhoods Af(x) have only internal grid-knots, Af(x) C int(G n (l, ±)). 
The associated diagonal submatrix (A n )jj is symmetric. Diagonal submatrices for 
which Af(x) contain boundary grid-knots of G n (l,±) may lack the symmetry. If 
the quantities a\j in ([36]) are replaced with ajjj , the symmetry of A n is lost 

altogether, although the convergence is still preserved. However, the quantities ffly 
must not be replaced with ciij(hk) since the resulting (A n ) k i would be discretizations 

of-YnjOiAdj- 



Figure 3: Numerical neighbourhoods at internal grid-knots of Di 




ai2 > ai2 < 



ri = 3, r 2 = 1 



4.2 Discretization with compartmental structure 

Now we can describe general structure of constructed discretizations of Aq(x). From 
the definition of bilinear forms a n (u,v) = (v|A„u) the following property is obvious: 
If v = 1 and u are with compact supports on G n there must be a n (l, u) = a„(u, 1) = 
0, implying that the row sums and column sums of A n have zero values. Hence, if the 
off-diagonal entries of A n are non-positive then the matrices A n are simultaneously M- 
matrices and have the compartmental structure. We consider here only the extended 
scheme. 

It is convenient to utilize the quantities: 

1 d 

— ^2 au{hk + hmii(l,s)), (38) 

where ma(l, s) are defined by the rules of construction of extended schemes. 

Discretization procedure 4.1 Let Assumption ^ ■ 1\ be valid and matrices A n on G n 
be constructed by the rule 131]) . Then their entries have the following properties: 

1. Entries of(A n )ki,k,l G I n , x = hk are linear combinations of a.ij(xij(n,x,l)) 
where Xij(n,x,l) — hk + hrrtij (I , s) , rriij(l,s) are n-independent elements o/ffi. d 
and I £ L, s G {1, 2, . . . , d}. 

2. For each grid-knot x — hk: (A n )kk = — ^2i(A n )ki- 

3. For each x = hk G cls(Di) entries on the "cross branches" x ± hei, i.e. 
{A n )k±h ei , are defined by: 

( A n) xx ± hei = -^[w»(o«,aO 

^ ^ 777" | ^im ip^im {p*i 0) I • 
TrnKf) 

4- Entries of A n which are not on the "cross branches" are defined by using ele- 
ments Zij(l) = ri(l)ei — rj(l)ej G /„ or elements Zij(l) = ?-j(Z)ej + rj(l)ej G I n : 

( A »)kk± ZlJ (i) = - n (l) rj (l) l a *i( aB «( n ' a! '0)l- 

Some peculiar features regarding the structure of sets N(x), x = hk G G n fl Di, 
must be pointed out. If a.y , i j is not changing the sign in a neighbourhood of x 
then the minimal number of elements in N(x) is 1 + d + d 2 . In this case the set Af(x) 
consists of its center, 2d-grid-knots on the d-dimensional cross {±e^ : i = 1, 2, . . . , d} 
and 2 grid-knots in each two-dimensional plane. Generally, the number of grid-knots 
in two-dimensional plane may be larger than 2. Here we consider only the case of two 
grid-knots at most in the two-dimensional planes. This demand has the following im- 
plication on the construction of discretizations A„ s \ Let the pairs e r ,e s and e s ,e t 



define two-dimensional planes and let An , A„ be the corresponding discretiza- 
tions which are constructed using the parameters r^ rs \r^ st \ Then there must hold 
(r( rs )) s = (r( st )) s . In such case the off-diagonal entries of A n have the structure as 
described in 4. of Discretization procedure 14. 11 The described structure of sets Af(x) 
is valid for all x because the functions dy , i ^ j do not change sign on S n (p, x) . This 
is an important consequence of the strict separation of sets max{ai2,0}. 

Obviously that all the constructed forms a n of this section are discretizations of 
the form ([3]). One can be easily convinced that the terminology "a discretization of 
the original form ([3])" is not artificial. At the present level of analysis it is easy to 
check a(v, u) = lim„ h d a n (v, u) for any pair v, u 6 

THEOREM 4.1 Let Assumption ^. 1\ be valid. There exist discretizations A n which 
are constructed by the rules of Discretization procedure \4-l\ such that A n have the 
compartmental structure. 

Proof: For each Di we have to choose the parameters r{(l) of the properties 
3. and 4. of Discretization procedure 14.11 so that the condition 3. of Assumption 
14.11 is valid. The uniform continuity of coefficients on Di and the condition 3. of 
Assumption 14.11 ensure the compartmental structure of matrices A n as demonstrated 
in [LR3] , QED 



5 Convergence in W^-spaces 

Discretizations of the original variational problem |5]) or (j6|) are defined in terms 
of a sequence of bilinear forms a n (v) on E n (R,M. ) x E n (R, R d ) and a sequence 
of linear functionals )jj on E n (R,M. d ). The associated discretized variational 

problems are defined by equalities (f27|> . The discretized variational equalities (|27|) 
can be rewritten in an equivalent manner: 

(XI + A n )u n = fj, n , (39) 

where in the case of Problem ((6|) u„, fi n are grid-functions on G n and in the case of 
Problem ([5|) they are grid- functions on G n (R, D) or G n (R, D). 

Problems ([5]) or © are solved numerically in two steps. In the first step we con- 
struct grid- functions u„ on G n (R) or G n (R, D) according to ([39]) . The obtained grid- 
functions represent the solution at grid-knots, therefore, we call them grid- solutions. 
The grid-solutions are imbedded into the spaces E n (R, M. d ) or E n (R,D) by (Tl6|) . and 
the convergence u(n) — > u must be proved in some Banach spaces. Therefore we call 
functions u(n) approximate solutions. 

Though the functions u(n) = $(i?)u n are called approximate solutions, this ter- 
minology has to be justified after a convergence analysis. The convergence proofs are 
based on some properties of the forms a„(-, •) and linear functionals ( - \n n )r to be 
described in details later in this section. Most of the analysis in this section is related 
to Problem ((6|) on M d . The obtained results can be easily applied to Problem ([5| on 
a bounded domain. This is carried out at the end of section. 



5.1 Consistency- 
Certain number of notions important for the convergence of approximate solutions is 
formulated in terms of sequences of functions with a particular structure: 

2J = {v(n) : neN}cU n E n (R,R d ), 
H = {u(n) : iieN}cUA(fl,4 



DEFINITION 5.1 (Consistency) The forms a n (; ■) on E n (R,R d ) x E n (R,M. d ) 
are consistent with the form if 

a(v,u) = lim h d a n (v(n), u(n)) 

n 



is valid for any pair 03, il of such that 03 converges weakly in W^QR ) to v, and 
il converges strongly in W2 (R ) to u. 

PROPOSITION 5.1 Let a sequence of matrices {A n : n G N} be constructed by 
basic schemes or extended schemes and let {a n (-,-) : n G N} be the corresponding 
sequence of discretizations of f3|). Then the forms a n are consistent with the form 

W- 

This important result is proved by a lemma which is formulated bellow. Let 03 and 
il be defined by (|40|) and converge weakly and strongly in W 2 (R d ) to v,u 
respectively. Expressions u> n (a,ij,x) for i = j are defined by (|38| . In this proof we 
extend this definition to the case i ^ j and define uj n (aij, x) — Oy (ajy (ti, x, I)), where 
Xij(n,x,l) = x + h m,ij(l,n,s) are constructed by the rules of basic and extended 
schemes. For each I G L and each pair i,j G {1, 2, . . . , d} there holds: 



lim / aij(x) (div(n))(x) (dju(n))(x) dx 

= lim uj n (a t j,x) I (div(n))(x) (dju(n))(x) dx, 

" x £ G„(P)nft( £ ) Jc n ( P (i), x ) 

Let us point out that the written identity is unchanged if we replace to n (aij, x) with 
a,ij(x). Let us consider a bilinear form on W 2 (M. d ) x IT 2 1 (R d ) defined by 

8 n (l,v,u) = V p(u>(ra,a:)) / d l v{y) dju(y)dy, (41) 

7^„„ JC„(r.x) 



x£G n (R)nDi 



where p is a uniformly continuous function on Di and •) is a transformation of 
R such that w{n,x) G C„(r,x). Then 

lim e n (l,v(n),u(n)) = / p{y) d l v{y) d j u{y)dy (42) 

for any pair 53, il of (f40|) . converging in W^R ) weakly to v and strongly to u, 
respectively. 

The object of next analysis is the bilinear functional on E n (R,M. d ) x E n (R, R d ) 
defined by: 

j n (l,v,u) = h d vo\{R) ^ p(w(n,x)) 5i(rih)v(hk) Bj(rjh)u(hk). (43) 
In particular: 

-f n (l,v(n),u(n)) = h d vo\(R) 2J p(w(n,x))v k iU kj , 

xeG n (R)r\D l 

where uj„j = Bi(rih)v(hk),Ukj = Uj(rjh)u(hk). 

LEMMA 5.1 Let £/ie sequence 03 converge weakly in W^f^L ) ^° w anc ^ ^ converge 
strongly to u in W}(R d ). Then 

lim j n (l,v(n),u(n)) = lim 8 n (l,v(n),u(n)) = / p(x) div(x) dju(x) dx. 
n " Jo, 

Proof: We have to analyze 6 n (l,v(n),u(n)) as n — > oo. To simplify notation 
we assume i,j G {1,2}. The indices A; = (ki, k 2 , . . . , k s ) are denoted shortly as 
fc = (fci,fe') and A; = (ki, k 2 , k"), k\,k 2 G Z. The corresponding r\,r 2 are denoted 
by r,s, respectively. After inserting expressions for v(R,n),u(R,n) into (|4Tj) and 
carrying out a straightforward calculation we get expressions: 

n (l,v(n),u(n)) = ^2 p(w(n,x)) p n {hk,v(n),u(n)), 

x£G n (R)nD, 



where p n (hk,v(n),u(n)) denotes an integral which for the case of i = j = 1 has the 
form: 



p n {hk, v(n), u(n)) 



/ d 1 v(R,n)(y)d 1 u(R,n)(y)dy = V (Vv|#) 

JC n (r,hk) k , v 



^2 Vkk'Uw (di1p k \dl^l) + ^2 Vk+rk'Ull' (dxipk+r\dx^i) 
l—k,k-\-r l—k,k-\-r 

= rhJ2(^k'Hi')(ai(rh)v)(hk)(a 1 (rh)u)(hk). 



k'l' 

Now we use quantities s k i — (V'felV'i)IIV'fcllr 1 an< ^ bring into mind their properties 
Ski > 0, J2i s ki = 1- 

9 n (l,v(n),u(n)) = h d vo\(R) ^2 P( w ( n , x )) s k'l' «(fcfe')i «(fci')i- 

Analogously we get for i — l,j = 2: 
6 n (l,v(n),u(n)) = h d ~ 2 vol(R") s k"i" p(w(n,hk)) J dz\dz 2 

k,l,k" ,l"£l n (R),x£Di J(k,r)x,J{l,s) 
V{k,l,k")i U (k,Ll")j 1pk(zi)lpt(Z2) + V{k,l,k")i U(k+r,l,l")j tpk+r(zi)^l(z 2 ) 
+ V{ k ,l+s,k")i u (kd,l")j V'fcO^lM+sO^) + «(fc,i+ s ,fe")i u (k+rd,l")j 1pk+r{zi)^l+ s (z 2 ) 

where J(fc, r) = [hk, h(k + r)] and i?" stands for the parameter set {r , r 3 , r4, . . . , r^}. 
Upon integration over z\ , Z2 we get 



1 



9 n (l,v(n),u(n)) = -h"' 2 vo\{R") 



^] s k »i»p(w(n,hk)) 

k,l,k" ,l"El n (R),iBEDi 

V(k,l.k")i U(k,l,l'')j + v (k,l,k")iU(k+rd,l")j + v (k,l+s : k")iU(k,l,l")j + V(k,l+s,k")iU[k+r,l,l")j 



We finish the proof for the case i = j = 1 since the case i ^ j can be treated 
analogously. The quantity j n would be equal to 6 n if s^i were absent and the double 
sum were replaced with the single sum over indices k. Therefore wc arc due to 
estimate their difference: 



j n (l,v(n),u(n)) - 9 n (l,v(n),u(n)) = h d vol(R) J~]p(w(n, hk)) v ki S k i s k >i> (uij-u k j). 

kl 

(44) 

To estimate the right hand side we need p = supp: 

7 n (l,v(n),u(n)) - 9 n (l,v(n),u(n)) < 
V h d/2 1 Ui {n)v n \ R2 h d ' 2 max{ | (Zfo , j) - I) U 3 (r,)u n \ R2 : j = 1, 2, . . . , d}. 

By Theorem 13 . 1 1 we have 

^It/^KIL < (1 -a 2 )- 1 \\div(R,n)\\l 

h d \{Z(r 3 ,j) - I) ^(r,)u„| 2 fl2 < (1 - a 2 )- 1 ||(Z(«i) - I) d j u(R,n)\\ 2 , 
where w = hrjej. Due to the strong convergence of il we have 



h<V*\(Z{r jt j)-I) 

< (1 -a 2 )- 1 ' 2 



U {r )u n \ m < (l-a 2 )- 1 / 2 \\(Z(w)-l)d J u(n,R)\\ 2 
\\{Z(w)-I)d jU \\ 2 + 2\\djU- djuiRrfh] -» 0, 



so that lim n j n (l, v(n), u(n)) — lhn n d n (l,v(n),u(n)) = J pdivdju. QED 
Now a proof of Proposition 15.11 follows from the inequalities 



lim V j n (l,v(n),u(n)) = / 

n — /i 

lee JU 



in which p(w(n, hk)) of P3"| is replaced with (n, x, I). 



(x) div(x)dju{x) dx. 



5.2 Strict ellipticity of discretized forms 

A discrete form a n (-, •) on lo(G n (R)) x l (G n (R)) is said to be strictly elliptic |Yoj if 
there exist two positive numbers M(a n ), M(a n ) such that 

d d 

Mtfln) E I 

Ui(n)u\ 2 R2 < a n {u,u) < M(a n ) ^ [Ui(ri)u\ 2 R2 . 

i=l i=l 

For a sequence of discrete forms a„(-, •) we need a stronger result. The strict ellipticity 
must be uniform with respect to n and vol(i?) = Y[ r i different values of Tq in the 
parameter set R = (r n ,r). 

DEFINITION 5.2 Let A n be discretized by the rules of Discretization procedure 
\4-l\ Discrete forms a n {-, •) on Zo(G„(i?)) x lv{G n {R)) are said to be strictly elliptic 
uniformly with respect to n £ N if there exist positive numbers M_ < M such that 

d d 

MJ2 I^KHL ^ a n {u n ,u n ) R < MJ2 l^(^)u|L- (45) 

i=l i=l 

for all n G N and all < (r )i < r,, i = 1, 2, . . . , d. 

PROPOSITION 5.2 Let the discretizations A n of A (x) = — ^ c^dy (x)dj be con- 
structed by the rules of Discretization procedure \4-l\ If A n have the compartmental 
structure then the discrete forms v, u t— > ( v|v4„u are strictly elliptic on lo(G n (R)) x 
lo{G n (R)) uniformly with respect to n e N. 

Proof: First we consider the basic scheme for a two-dimensional grid. From the 
structure of bilinear forms ([32]) we have 

a n (u,u) > MAUxipdunW + M(\K{-)U 2 {-P2)u n \l+ \K(+)A 2 U 2 (-p 2 )u n \fj, 

(46) 

where K ( =p ) = 11 g „ ( p, t ) are pro j ectors . Let us consider a pair x , y = x — hp 2 e 2 which 
is involved in the definition of operator U 2 {— p 2 )- The indices of x, y are k, I, respec- 
tively. If x,y £ G„(P, +) then {U 2 {—p 2 )\i n )k can be replaced with {U 2 {p 2 )u n )i. If 
y £ G n (P, +) then y s G n (P, — ). In this case the term (U 2 (—p 2 )u n )i can be omitted 
from the sum in (|46[) because this term is already contained in the corresponding sum 
of Ui(pi)u n . Hence, the right hand side of (|4"6"| can be estimated from bellow by 
JkLY^=i \Ui(pi)u n \2- In this way the left hand side of inequality (|45|) is proved. The 
right hand side of this double inequality follows by choosing the double value of M. 
In the case of d > 2 we use the construction (pTTj) and get the same lower and upper 
bounds. Let us point out that the compartmental structure is not used in this step 
of proof. 

Let us now consider a two-dimensional problem with an extended scheme. In the 
present step, the compartmental structure is utilized in an essential way. The form 
a n ~\u,u) of (|34|) can be rewritten as: 

a„~\u, u) = 

E E E (^(^ + < (++) (l)) + ^«i 2 ( a; + t (++) (r)))(a i (/ l )u)( : r) 2 - 

J2 E «i2^+* (++) (r)) (^(ai^^w-^^t^ww) ■ 

(47) 

Let us remind that a\ 2 < on the set D{—). Due to the compartmental structure 
the first term is positive definite. The second term is positive semidefinite and can 
be disregarded in the next step of estimation from bellow. The result is: 

al\u,u) > u{a) E E E ( U *Ol 
leC- keG n (i~) i 



where uj(a) is the positive number specified in Assumption ^. 11 For the form a„ ( u j u) 
we have an analogous inequality involving the summation over all the indices k £ 
G„(+). 

al +) (u,w) = 

E E E(^^ + t(+_) ( 1 ))-^ a i2(a ; + t( + -)(r)))(a i ((-ir 1 / i H( : r) 2 + 



E a 12 { X + t^-\r))( x h^^ 1 {r 1 h)u)(x)- 
iec + keG„(i,+) \V r2 ^ 



" :(0 -(a 2 (-r 2 / l ) u )( a; ) 



n(0 

(48) 

The lower bound follows in the same way as foral _) . QED 
5.3 VK^-convergence 

We have shown how the form A(v|u) + a(v, u) is discretized by forms h d X( v n |u n }r + 
( v n \A n u n )r. In order to solve discretized problem ([55)1 , we have to describe a dis- 
cretization of the linear function v —> (v\n) by h d ( w n \fji n )r, where fi n £ l(G n (R)). 
First we must demonstrate the existence of n n such that h d ( v„|/x„ )r — > 
Discretizations of /x are defined by ([21]) so that is valid. 

LEMMA 5.2 Let /i oe a continuous linear functional on WqCR )• There exists dis- 
cretizations fi n (R) £ l(G n (R)) such that 

(u(n)\fx) = /i d (u n |/x„) fl = /i d vol(i?) u ^k- (49) 

kei n {R) 

Obviously, for each sequence {u(n) : n £ N} weakly converging to some u £ W% (R ) 
the following equality holds: lim n (u(n)|^t) = 
From this Lemma we have 

l(u„|/Ofll < ft-^VK) 172 ll^l| a ,-i. 

Inequalities ([43]) . (j4T)|) and the variational equalities (|27p imply the first result 
towards our proof of convergence of approximate solutions. If u n solve (|27|l or 
then 

^MKHj < /i d (u n |(A/ + A„)u n > fl < HJe,n)|| 3 ,i||/i||2,-i. (50) 

COROLLARY 5.1 Let u n = T(A, and u(i?,n) = $„(-R)u„. T/ien /or eac/i 

i? i/ie sequence il = {u(i?,n) : n G N} C U n E n (R, R ) converges weakly in W^R^) 
to some u G W^(R d ). 

Let it* be the solution of (J6j> . Then the sequence of functions u*(n), defined 
by (j20|) . strongly converges to u* in W2 (R )• In the remaining part of this analysis 
we have to demonstrate the expected property lim„ u(R, n) — lim„ u*(R, n) — u* for 
each R. We follow the well-known finite element technique. 

Mh d \u n -K\m,i < h d (u n -u* n \(XI + A n )(u n -u; i ))R 

= h d (u„ - u* I (XI + A n )u n ) R , s 

- /i d ( u„ - u* I (A/ + A„) u* } R [ ° > 

= h d ( u„ - u* I fi n } R - h d (u n -u* n \ (XI + A n ) u* ) R . 

By Lemma 15.21 the first term on the right hand side converges to (u — u*\fi). By the 
consistency property of Proposition 15.11 the second term converges to the same value. 



THEOREM 5.1 Let 51 be as in Corollary I5.il Then the sequence 11 converges 
W\ (M.)- strongly to the unique solution u* to (G|). 



From this result, Lemma 15.21 and Lemma 13. II we get another important result for 
A = 0. 



COROLLARY 5.2 Let D be a bounded domain with Lipsithz boundary and fi 6 
W^' 1 (D). Let A n (D) be the restriction to G n (R,D) of A n , /x n on G n (R,D) satisfy 
and u„ = A„(D) -1 ^i ra . Then the sequence it converges strongly in W%(D) to the 
unique weak solution u of {5|j. 



6 Convergence in Li-spaces 

In this section we consider Problem |5J for a bounded domain D and numerical so- 
lutions in Wp(D). Furthermore, for the sake of simplicity of exposition, we restrict 
our analysis to the differential operator Aq{x) = — Y)^ didij(x)dj and its discretiza- 
tions A n (D). We consider the boundary value problem (O with p e Tl{D) and its 
discretizations 

A n (D) u„ = n n (52) 

on G n (R, D), where A n {D) have the compartmental structure. The discretizations 
fj, n £ l(G n (R,D)) are defined by 

Pk = n , ii [Vk \p)- 

Apparently, (v(n)\fi) = v„|/x„ for any v(n) — J2 v k' t l J k S E n (R,D). Problem 
([52) is defined for grid-functions on G n (R,D). 



6.1 Boundedness of approximate solutions 

For any pair r, R, < r < R < 1 and v 6 D we define the balls B r (v), Br(v). By us- 
ing the functions it £ W^iD) we define the sets A(r,s,v) = {x 6 B r (v) : u(x) > s}, 
where sgR. The measure of A(r, s, v) is denoted by a(r, s, v) = meas (A(r, s, v)). Let 
us assume that there exist a subset S C W^C-D) and two numbers, c\, C2, independent 
of r, i?, s, t), such that the inequality 

d 

^2\\^A(r,s,v)d t u\\l < Ci || t A (R,s,v) fas) ||1 

2 = 1 

(53) 

is valid for a fixed q > d, each u S S, and all s > so with some sq. Then |LU| the 
functions it G S are bounded on D with a bound which depends on D, ci, C2 and ||u||2- 
Instead of balls B r (v) one can use rectangles S n (p, v) = Yli=i[— hpi + hv i} hvi + hpi] 
as well. Actually, it is sufficient to consider balls (rectangles) with centers v in a 
dense set of D whose radii (edges) are contained in a sequence {r m : m G N} C (0, 1), 
converging to zero. For instance, the sets S n (p, v) and v g G„(i?), match this weaker 
condition. This fact enables a straightforward application of ([53]) to the sequence of 
functions in Corollary 15.21 First, we have to establish a discretized version of (153]) . 
and then we have to prove that the constructed discretized version implies (|53|) for 
the sequence il of Corollarv l5.2l Further, let 

A(r,s,v) = S m (p,v) n {u(x) > s}, 
A(R,s,v) = S m - t {p,v) n {u(x) > s}, 

where the sets A(R,s,v) are larger than A(r,s,v), dist(<9A(r, s, v), dA(R, s, v)) = 
h(m)(2 t — l)p, where p — minpk and m > t. Hence, the symbols r, R stand for the 
d-dimensional parameters h(m)p and h(m — t)p 1 respectively. The discretization of 
A(r, s, v) is defined by F n (r, s, v) — G n (R, A{r, s, v j). The index set of F n (r, s, v) is 
denoted by J n (r, s, v) and its cardinal number by j n (r, s, v) = card( J„(r, s, v)). 

Let r < R e N and let us define the cut-off function 8(r, R, •) on R as a con- 
tinuous piecewise linear function, such that 6(r,R,x) — 1 for x G [— r,r] , and 
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zero outside of [-R,R]. Apparently, 6(r,R,-)' = (1/(R - r))l { _ R _ r] - (1/(R - 
r))lr ri ffl. Let r = (ri,r2, . . . ,r*d) and = (i?i,i?2, . . . , ■ The continuous function 
x i — ► x(r, i2, a;) = fL = i ^C 7 *) R»> l^i — v i\) ^ s a cut-off function on R d for which 
|5x(r,H,u, -)| < maxi(i?i - n) -1 . 

Discretizations of x(f", -R, u, •) on G n (R) are denoted by Xn( r > R> ' u ) an d they are 
defined in terms of rougher grids G m — t C G m C G„ for which m — t < m < n. For 
each 17 e G n (R), we define three sets 

if(r) = S^Op, v) = {a; : x(r, Jl, v, x) = 1}, 

tf(Jl) = S m _ t (p, «)=supp(x(r,H,t>)), (54) 
= K(R) + S n (p,0) = U veK{R) S n (p,v), 

which are illustrated in FigureHJ The rectangles (|54"|) have edges 2r,j = p,j2 1_m , 2i?j = 
Pi2 1+t ~ m and p i (2 1+t_m + 2 1 ~") units, respectively. It is important to notice that 
the grid-functions X„(r, R, v) arc defined by sets S n (p,v) which are related to the 
numerical neighbourhoods of constructed schemes. In our next proof we again use 
the fact that the sets supp(max{ay , 0}),i ^ j are strictly separated, so that we can 
assume that the functions ay on rectangles (|53jl do not change sign. 

It is easy to verify that \Ui(pi)x n ( r : R, v)\ < (h(m)p(2 t — l)) -1 , so that we can 
write 

\U t ( Pl )x n (r,R,v)\ < 1 h{n? 1 - - 2 < 

where L = 2 t+n ~ rn 1 l = 2 n ~ m , and p is a number depending on p = minpj(Z). The 
sets F n (r, s, v) and F n (R, s, v) are defined in terms of the sets (|54|) . 

F n (r,s,v) = K(r) nsupp(w„), 
F n (R,s,v) = K + (R) n supp (w„). 

The following estimates are used in our next proof: 

max Z(i,± Pi (l))x(r,R,v) < 1f(r,s,v), 

P (55) 
max \Ui(pi(l)) (Z(i,±pi(l))x n (r,R,v))\ < _ 1 F (r, s ,v)- 

Now we consider the following auxiliary problem: 

d 

A n (D)u n = f Q „ - ^2 Vi(pi)fin ■= A*„, (56) 

i=l 

where u„ are grid-solutions and U are grid- functions on G n {R, D). 



LEMMA 6.1 Let u„ = A n (D)~ 1 fi n , where \x n are defined by i56\) . There exist 
positive numbers Ci,C2 and Sq, independent ofn,v such that the inequalities 



\x n (r,R,v)U t ( Pl )u\ 2 R2 < c x \\ F{ R,s,v)( u n s)\ 



2 

R2 



are valid for s > sq and q > 2. 



C2S 2 J2 \^F(R, S: v)tin\R q (j(R,S,v)) /<? 

(57) 



Proof: Because of the rule (|3Tj) of construction of discretizations for d > 2, it is 
sufficient to prove (|57|) for two-dimensional problems. Basic schemes are considered 
first. 

The function w — max{u — s. 0} has the discretizations w„ and w n £/iW„ = 
w„[/iU n . Now we have to evaluate the forms 

a { n \x 2 w,u) = ( x n (r, R, v) 2 w n \ A [ n ] u„ ) R , ^ 
a { ^\x 2 w,u) = ( x n (r, R, v) 2 w n \ A^' u n } R , 

in terms of matrices Jl n ■ (i, j, n, rj) of Expression ([33]). It is sufficient to consider the 
form a n and one of the terms. For instance, the term involving the matrix An ' (I) 

is ~Ey=i K(Pt)^n +) (i ) i,Pi,Pj)U 3 -(Pj)- 

In order to write expressions in this proof in a concise form we omit various indices 

in the notation. Thus we use u,w, x,Ui and A^\i,j). Apart from this, we need 

the notations Zif = Z n (pi, i)t n , so that we can write Uif n — {hp{)~ 1 (Zif — f). In the 

equality C/j(gw) = (Uig)Zi\v + gUiw we insert g = X 2 to get: 

CA(x 2 w) = X 2 + (U iX ) (Z iX + x) ZiW. (59) 

Also we need 

{Z i x + x)U j w = 2 X U j w + ^(U k x)Z j w-^(U k x)w. (60) 

Pj Pj 

In this way the finite differences UjW are always multiplied by the grid-function x 
(not with Zkx)- 

We start the proof by an analysis of (|58p : 

<X 2 w|4-)(Z)u) R = J2( xU i v f \A(+\i,j) X U j w) R 
+ ( (^ w ) UiX\A (+ \i,j) (Z lX + x) U j w) R . 

The first term on the right hand side can be estimated from bellow as a consequence 
of the strict ellipticity (jl5|) of Proposition 15.21 The left hand side of this equality is 
equal to (x 2 w, \fi) R) so that we can write: 

1=1 (61) 
- J2( Z i wU iX\ A(+) (hj){ZiX + x)U j w) R . 

The two terms on the right hand side are denoted by Ti, i = 1, 2, respectively. Let us 
estimate them from above. The representation of /x, in (|56[) . is inserted into the first 
term to get: 

2 

|Ti| = \(x 2 ™\v)r\ < |<X 2 w|fo>fl| + K^(x 2 w)|f 4 ) fl .|. 



First (f59|) and then (|60[) are utilized in the last two terms on right hand side. The 
functions x, Zj~X are estimated from above by (I55|) and then the CSB inequality and 
the inequality ab < (e/2)a 2 + (l/2e)b 2 with a convenient choice of e, is applied to all 
the terms on the right hand side. 

2 

\T X \ < e]T \ X Uiw\%2 + + fR— jqa) I^CH.-.*) 

^ ' i=0 

where &i does not depend on n, e. It is important to notice that the first term 
on the right hand side of this inequality and the left hand side of (IrJTj) involve the 
same norms, and can therefore be subtracted. Now we use the Holder inequality 
|iF(R,s,«)fi||',2 _• \^F(R.s,v)(i\ji q j(R, s^) 1 " 2 / 11 and get the first intermediate result: 

2 

\Tx\ < eJ2 \xUiw[ 2 R2 + h (l + -ja) |lF ( ^)w| 2 fl2 

i=l 1 1 

2 

+ 62 E |lF(H,-,«) f i|H,i(«»'»»*) 1 " 2/ *. 

where 62 depends on 1/e. 

The quantity T 2 can be estimated analogously. Let us use (|50j) in Expression (|6"Tj) 
and apply the same technique as in the previous case. We get: 

2 1 
\T 2 \ < U^w| 2 fi2 + b 3 |R _, 2 |lf(fl, a ,«)w|^, 

4=1 ' ' 

where 63 depends on M, 1/e. Thus the inequality (|61j) implies the following result: 
2 

(M-2e) J] IxMSq < h (l + | fl _ r | 2 ) |lf(H,.,«Hffl 

i=l ' ' 

2 

+ &2^ |ljr ( «, a ,„J ft |St, J(JI, *, w) 1_2/ «, 

8=0 

which is equivalent to the assertion of lemma. 

Now we consider the extended schemes. Let the discretizations A n (D) of |57|) be 
constructed by extended schemes and let u„ be the corresponding solution. Then we 
have a n (v,u) — (v n |/i„), where a n (v,u) is defined in (|34[) . Let us define v„ = x 2 w 
and let its image v(n) = $„v„ be inserted into the expression a n (v,u). If we apply 
(|59|) to the function v(n) we get: 

a n (v(n),u) = a n (x;w,w) + ^6 n (s,x,w), (62) 

s 

where the terms on the right hand side are defined as follows. The form a„(x; w, w) is 
the sum of forms (|T7|) and in which the factors (a; + (r)) are replaced with 
the factors x( x ) 2a ij( x + * E acn &n(s,X,w) is a form which is a first degree 

polynomial in Ui(pi(l)))w ni with various l,i,pi(l), resulting from the application of the 
rule (l59l) . As in the corresponding proof of Proposition 15. 2, the quantity a n (x;w,w) 
can be estimated from bellow, 

2 

a„(x;w,w) > u(a) |xC^w|l. 

i=l 



An attempt to estimate the remaining terms on the right hand side in (f6"2")l from 
above as in the first part of proof would lead to an unresolvable problem. Among 



the resulting terms there would be sJ2i |%«^i w n|i as wen as £ J2i \Xn^i(Pi) w n !!• 
The former one can be moved to the left hand side of (j62|) and subtracted from 
a; ( a ) Si |Xn^i w ra|l- Unfortunately, this cannot be done with the latter one. There- 
fore, before finding upper bounds on & n (s,x,w) we have to apply another version 
of (5Q): 

(x„ + Z„(pi, i)Xn) u j(Pj)™n = 2x n U j (p j )w n + 

Pi 1 

— (Ui(pi)Xn){ Z n{Pj,j) - J)w„ H (7 + ^ n (ft,i))(x„ - Xn) ( Z n(Pj,j) ~ ^)w„. 

Pj L J 

If these expressions are applied, then estimates from above of b n (s, Xj w ) contain the 
terme^ |x n ^i(pi)w„|f. Let us choose x(x) = x(r, R, v, x)x(r, R, v, x+h(p%ei + 2 
e 2 )). Then 

su PP (x) + C„(p,0) = K(R), 

and (x n C^(Pj) w n)fe is a linear combination of (f7j(l)w n )j such that hi £ K(R). Now 
we have 

s ^2[Xn U i(Pi) w n[% < E X! |Xn t/ *( 1 ) W n I2 

and the difficulty regarding |x n £/i(pi)w n || with various ^ 1 is removed. The 
grid- functions h-\(Xn ~ Xn) nave estimates from above as in (|55|) . In this way we 
get (Jill) again. QED 
The function w(n) — max{u(n) — s,0} and the grid-function w„ which is compo- 
nentwisely by Wf, — maxjufc — s, 0}, are related by the expression w(n) — $„(i?)w„. 
The index set of grid- knots of G n (R) in the (closed) set S m (p,v) is denoted by 
I(m,p,v). The restriction of w(n) to the set S m (p,v) and the function w(n,I) = 
Efee/(m p v) w k'4'k ar e not equal. They are equal on the set S m (p, v) and the following 
rough estimate is valid elsewhere: 

Mn,I)\\t = \\l Sm{PtV) w(n,I)\\l + \\(l-l Sm{p , v) )w(n,I)g 
< 2 d \\l Sm{p ^w(n,I)\\l = 2 d \\l Sm(p ^w(n)\\l 

Now we combine the result of Theorem 13. 11 i.e.: 

(! - o*) \li( m ,p,v) w„|| 2 < IIM^-Olll ^ |il(TO,p,«) w„l| 2 

and the previous inequality in order to get: 

1 2 d 
[li(m,p,v)Wn\ R 2 < 1(j2 \\w(n,I)\\l < — \\1s m { P ,v)w{ri)\\ 2 (63) 

It remains to compare the functions diw(n) and E/j(p,)w n . By Expression (|26[) we 
easily prove: 



Us m (p,v)9iw(n)\\l < || T,kei(m. P ,v) diw(n)\\l < h d |l/ (m , p ,„) (U l (p l )w n ) J 2 R2 
< ^|x n (r,H,«)fi(Pi)w n |f 



I 2 



(64) 



(65) 



Let us multiply ([57)1 by use and afterwards The result is: 

2^ ci / 1 

Ef=l UA(r,s,v)diw(ri)\\% < 1 _ ^ 2 ||lU(H,s,„) ( 1 + | fi _ r | 

+ c 2 s 2 (h 2 ^J2 \tin\%,)a{R,s,v) 1 -*'*. 

i=0 

In this way we obtain an important property of grid-solutions to Problems ([5; 

COROLLARY 6.1 Let D C R'' fee a bounded domain with a Lipshitz boundary and 
q> d. If h d ^2 d =0 | fin||j g are bounded by a number uniformly with respect to n, then 
the grid- solutions u n on G n (R,D) to A56}) are also bounded uniformly with respect to 
n. 

Proof: One has to compare (p)3"|) and inequality (|65|) . QED 



6.2 Weak consistency and convergence 

As in the case of W / 2 1 " a PP roacn i the convergence in L p -spaces is proved by utilizing 
certain kind of consistency. This weaker consistency is defined in terms of sequences 
similar to l[30]): 

SU(p,c) = {v(n) : neN, h d Y,i\UiV n \ p p < c} C U„£„(R d ), 
Ho = {u(n) : n G N, it £ C^ 3 (M d )} C U„£„(R d ). 

DEFINITION 6.1 (Weak consistency) We say that forms a n (-, •) on S n (M d ) x 
EV^R'') are weakly consistent with the form |3J) i/ 

a(v,u) = \im h d a n (v(n),u(n)) 

n 

is valid for any weakly convergent 2J(p, c),pG [1, oo), c > and any ilg o/ {5; 



LEMMA 6.2 Let a sequence of matrices {A n : n 6 N} be constructed by basic or 
extended schemes, and let {a n (-, •) : n G N} be the corresponding sequence of dis- 
cretizations of |3J). Then the forms a n are weakly consistent with the form |3J). 

Proof: The present proof and proof of Lemma HTTl are the same up to the equality 
(|44|) . The obtained equality must be estimated now by the Holder inequality: 



j n (l,v(n),u(n)) - 9 n (l,v(n),u(n)) 



< 



P h d ^\Ui(ri)v n \ Rp h d / q max{| (Z n (r Js j) - iJl/jfaKl** : j = 1,2, . . .,d}. 

Thus we have to consider | (Z n (r -, j) — l)Uj(rj)u n \ Rq for a fixed jeN and prove that 
this quantity converges to zero as n — » oo. It is easy to get the following expression: 

((Z n (r j ,j)-l)U j (r j )u n ) k = ■^^(i> h \{Z(the j )-l)B j (r j h)u) = (wn(r,-)) fe , 

where the function iv(rj) = (Z(thej) — i) 3j(rjh)u is continuous and converges point- 
wisely to zero as h — ► 0. By (ii) of Lemma we get 



h d/q \^n(r 3 )\ Rq < \\w( rj )\\ q — 0. 

h -> 



QED 



THEOREM 6.1 Lei £> C M d fee a bounded domain with a Lipshitz boundary, and 
let H = {u(n) : n G N} C E n (R, D) be a sequence of approximate solutions, u(n) — 
$, i (_R)u rl; where u n are grid-solutions to i52\). For each p G [l,d/(d— 1)) there exists 
a subset J{p) C N such that il = {u(n) : n G J{p)} converges to the unique solution 
to ([5p strongly in L p {D) and weakly in W^{D). 

Proof: The proof is split into two steps. In the first step we have to prove that for 
each p G (1, dj (d — 1)) there exits u G Wp (D) and a subset J(p) C N such that 

u = w— lim u(n) G W^(D), 
neJ(p) F 

u = s— lim u(n) G L p (D). 

In the second step of the proof we demonstrate that u coincides with the solution u* 
of the original problem ([5]) . 

Because of D C S\(t,v), for some t, v we have ||tt[[ p < (maxfcifc) (max/c ||9fctt||p) 
so that we can consider only d\u{n) and the corresponding U\(pi)u n . By (|26|) and 
estimating the function ip^i by the indicator of its support we get: 



\d lU (n)\ < 2 d \( U M 
hei n (R) 



Figure 5: Domains for Examples 17.11 and 17.21 
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where v = hk. Hence 

\\dMn)\\ p p < 2<^|£/ 1 (p 1 K|V 

It remains to demonstrate that the right hand side is bounded uniformly with respect 
to n. We can write 

|[/l(pi)u n |_R p = (V„ \ Ul(pi)VL n ) R> 

whevev k = / l -^ sign (^ 1 (p 1 )( Un ) fc )|^ 1 ( :pi )( Un ) fc |P-i/| C / 1 ( : p 1 ) Un |P-i. Then/i d /9|v„|^ = 
1. Therefore, 

|£/i(pi)u„| Rp = - (Vi(pi)v„ \u n ) R = - (A" 1 (Vi(pi)v n ) | fi n } R . 

By Corollary 16. II we get |A~ 1 T4(pi)v„ l^oo < (5 with some (3 depending on D and 
the bounds of h d | v„ | q Rq . The final result is h d /p\U l {p i )u n \ Rp < (3. 

The unique solution to (JSJ is denoted by u*. For p g (1, d/{d - 1)), v € 
and the corresponding v„,{>(n) = $ n (i?)v„, the following relations are valid: 

1) v = s-lim«(n) E Wl(D), l/p+l/q=l, 

2) /i d (v„ l^)^ = (£(n)|/z) -» {v\n), 

3) lima(w(n),u(n)) - lim h d (v n \A n (D)u n ) R = 0. 

n n 

The inclusion 1) follows from Lemma 13.41 and the uniform continuity of dv\ on D. 
Although this inclusion is plausible, some technical details are needed in order to 
transfer d\ from div{n) to (div){n). The relation 2) follows from v € C(D) and the 
construction of fx n . The identity 3) is the weak consistency which is proved in Lemma 
16.21 Notice that the pair v, u in this proof and the pair v, it of Lemma 16.21 have roles 
interchanged. Now we have 

a(v,u) — a(v,u) — a{v{n),u{n)) 

+ a(v{n),u{n)) - h d ( v„ | A n {D) u n ) R 
+ h d (v n \A n (D)u n ) R - h d (v n \v n ) R 
+ h d (v n \n n ) R . 

The right hand sides in the first three rows are either zero or converge to zero as 
n — > oo. Thus a(v,u) = (v\n). Because C«(L>) is dense in W^{D),p e (l,d/(d-l)) 
we get u = u*. However u - u* = <E WUD) implies u - u* = E Wl(D). QED 



7 Numerical examples in two dimensional case 

We test the numerical methods on examples for which we know solutions in the closed 
form. In the first example we solve the problem Q with an — and a measure on 
the right hand side, and in the second exmple we have a piecewise constant a\i and 
a non-trivial application of the extended scheme. 



Example 7.1 Let the matrix-valued function x i— » a(x) be defined on R 2 by 



a(x) 



1CT 4 x 



1 
1 
a 2 
1 



for xi < 0.5, 
for x\ > 0.5, 



and let us consider the elliptic differential operator A(x) = — Y^j=i di®ij(, x )dj on 



For t = (0.5, 0.5) the differential equation A(x)F(x 
F(x) = -c 



6(x 



In ((xi - 0.5) 2 + (x 2 - 0.5) 2 ) 
'ln((xx -0.5) 2 ) + a 2 (x 2 -0.5) 2 ) 



t) on R 
for 



ln(o- 2 ) 



for 



has a solution 

xi < 0.5, 
Xi > 0.5, 



where 



= 2tt 1 



The function F(-) is continuous on R 2 \ {£} and its first partial derivatives have a 
jump at x\ — 1/2. Let D = (0, l) 2 C R 2 as illustrated in Figure[5] In this example we 
consider the boundary value problem A(x)u(x) = S(x — t) with the nonhomogeneous 
boundary conditions u\dD = F\dD. Obviously, the solution u coincides with F on 
D. 

In numerical calculations we used a 1 = 10. Hence, we have two sets for the 
discretization procedure 14.11 D\ = {x € R 2 : x\ < 0} and D 2 = {x 6 R 2 : x 2 > 0}. 
The extended scheme was utilized for the discretization of A{x). Discretization of R 2 
is realized by grid-knots x = hk\e\ + hk 2 e2, h — 1/400, so that D is discretized by 
399 x 399 grid-knots. The system matrix A n of (|39| has the order 399 and the linear 
system is solved iteratively: if K n = diag(^4„) and Q n — diag(A„) — A n , then Q n > 
and 

oo 

{K- Y Q n ) k K-\ (67) 



A: 



E 

k=Q 



Let u n (r) be the approximation of u„ after r iterations. By taking the stopping 
criteria to be |u„(r + l) — u„(r)|i < 10~ 9 , we have found that the iteration terminates 
after r = 182000 iterations. Then we compared the numerical approximations u„(210) 
and values of the solution in l aD (G n (D) \ {t}) and l\(G n (D) \ {i})-norms, where the 
set G n (D) \{t} is chosen naturally instead of the set G n because the solution u of this 
example is not defined at t. Let the corresponding norms be denoted by | • |' ,p = 1, oo. 
The solutions are compared according to the relative error: 



|u* - u„ (210)|; 



|u*| 



p= l,oo, 



(68) 



where (u*) fc = u(hk). We obtained £i jre ; = 0.003371 and £oo,rei = 1-6254. The value 
£oo,rei = 1-6254 is realized at the grid- knot t—hex, i.e. at one of the nearest grid- knots 
to the singular point t of solution. Actually, the difference between two solutions at 
this grid-knot is 0.0075. Let us mention that the respective error E\ :re i = 0.267 is 
obtained with the grid-step h = 1/200. 



Example 7.2 Next, let us consider a differential operator A = — ^i a ijdj with 

the diffusion tensor 



a(x) 



a(x) 



t{x) = pt Do {x), p 2 < (T 2 , 



where a 1 is a positive number, p is a real number and Do = (1/4, 3/4) 2 (see Figure 
[S]). The function x t— » u*{x\, x 2 ) = x\X 2 is the unique solution to the boundary value 
problem 

(Au)(x) — p{x) for x E D, 
u\dD = u*\dD, 



where 

H(x) = 2pt Do (x) + f [<5(^i - |) -3*(xi-|) + *(*2-^) -3*(aa-|) 



The set K 2 is discretized as in the previous example. The sets G„(2) = Do n G„ 
and G„(l) =G n \ G„(0) define a partition of G„. Let G„(l, £>) = G„(l) n -D. Then 
G„(2), G„(l, £>) is a partition of G n (D) to be used in constructions of numerical grids 
and approximate solutions. 

To demonstrate the efficiency of the extended schemes we choose a 2 — 10, p = 2 
and the scheme parameters r\ = l,r2 = 3 as illustrated in Figure [3] These val- 
ues of parameters ensure the compartmental structure of the system matrix. The 
linear system (|52l) is solved iteratively as in the previous example. Let u„(r) be 
the approximation of u„ after r iterations. By taking the stopping criteria to be 
|u„(Y + l) — u„(7-) 1 1 < 10 -9 , we have found that the iteration terminates after r = 210 
iterations. Then we compared the numerical approximations u„(210) and values of 
the solution in i 00 (G n (.D)) and Zi(G ra (-D))-norms. We obtained ei yre i = 0.8 x 10~ 5 . 
In addition, |u* — u n (210) |oo,re; = 0.4 x 10~ 2 is realized at the grid-knot with coor- 
dinates x = (0.75,0.75). 



8 DISCUSSION AND CONCLUSION 

Attempts to discretize a second order elliptic differential operator by matrices with 
compartmental structure go back to the work of Motzkyn and Wasov [MWj . Their 
construction is based on rotations, which can be avoided in 2-dimensional cases by 
using extended schemes which we propose in this work. These can be used to get 
monotone schemes for any 2-dimensional problem with the second order elliptic oper- 
ator, in divergence or non-divergence form. For operators in divergence form one has 
to use the construction of system matrix described in Section [H while for operators 
in non-divergence form system matrices are of a simpler structure. Entries of system 
matrices are linear functions of dij(v), where v is the considered grid-knot. 

In the case of dimension d > 3 the construction of Section [4] does not always 
produce monotone schemes. An additional condition on the matrix-valued function 
x 1 ► Wij(x)}ff which is described in Assumption 14.11 assures the compartmental 
structure of matrices A n . 

LEMMA 8.1 Let A n on G„ be consistent discretizations of the differential operator 
A(x) — —^2ijdiaij(x)dj and let the entries [A n ) kl be linear combinations of aij(z) 
with some z £ S n (p(l), hx). Then A n can possess the compartmental structure iff the 
tensor-valued function a is strictly positive definite on M. d . 

Proof: For a constant diffusion tensor we prove easily this lemma by demanding 
the following property at each grid-knot: 

(A n i n ) k = A(x)f(x), x = hk, 

where x 1— > f(x) is any second degree polynomial in variables Xi. This property is 
equivalent to the consistency. To prove the necessary part, we assume that a(x ) is 
negative definite and A n have the compartmental structure. Then the same assertion 
must hold for a tensor- valued function which is constant locally at Xq . This contradicts 
the first step. QED 

Hence, if the matrix-valued function a is indefinite the only way out is to use 
rotations. Unfortunately undesired features may appear; in a subdomain Di C D 
where a rotation is required, some larger diagonal submatrices can be reducible. This 
causes the surfaces u(n, x) = const to have a saw-like behavior. In this case an 
averaging procedure of grid-solutions helps to obtain reasonable results. 

A different approach to the construction of monotone schemes for a class of elliptic 
operators was presented in Samarskii et al. [SMMMj . They constructed schemes with 
various grid-steps along coordinate axis which we called basic schemes in our analysis. 
In this way they obtained monotone schemes for a wide class of elliptic operators. 
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